Realization of the mean-field universality class in spin-crossover materials 
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In spin-crossover materials, the volume of a molecule changes depending on whether it is in 
the high-spin (HS) or low-spin (LS) state. This change causes distortion of the lattice. Elastic 
interactions among these distortions play an important role for the cooperative properties of spin- 
transition phenomena. We find that the critical behavior caused by this elastic interaction belongs to 
the mean-field universality class, in which the critical exponents for the spontaneous magnetization 
and the susceptibility are (5 = 1/2 and 7 = 1, respectively. Furthermore, the spin-spin correlation 
function is a constant at long distances, and it does not show an exponential decay in contrast to 
short-range models. The value of the correlation function at long distances shows different size- 
dependences: 0{1/N), 0(l/\/iV), and constant for temperatures above, at, and below the critical 
temperature, respectively. The model does not exhibit clusters, even near the critical point. We also 
found that cluster growth is suppressed in the present model and that there is no critical opalescence 
in the coexistence region. During the relaxation process from a metastable state at the end of a 
hysteresis loop, nucleation phenomena are not observed, and spatially uniform configurations are 
maintained during the change of the fraction of HS and LS. These characteristics of the mean-field 
model are expected to be found not only in spin-crossover materials, but also generally in systems 
where elastic distortion mediates the interaction among local states. 

PACS numbers: 75.30.Wx 75.50.Xx 75.60.-d 64.60.-i 



I. INTRODUCTION 



Spin-crossover (SC) materials consist of local units 
(molecules), each of which has two different spin states, 
i.e., the low-spin (LS) and high-spin (HS) states. The 
LS state is energetically favorable and dominates at low 
temperatures, while the HS state dominates at high tem- 
peratures because it is entropically favorable. The tran- 
sition between the LS and HS states is also induced by 
changes of the pressure, magnetic field, light-irradiation, 
etC'^'^'^i^i^'^'^ When interactions between molecules are 
weak, the HS fraction changes smoothly with tempera- 
ture. However, when the interactions become strong, the 
system exhibits cooperative phenomena.® The change in 
the HS fraction becomes sharper with increasing inter- 
action. When the strength of the interaction exceeds a 
critical value, the change becomes discontinuous. In or- 
der to control electronic and magnetic properties of SC 
compounds, it is important to understand the bistable 
nature of such molecular solids. 

As an important ingredient of the spin-crossover tran- 
sition, we need two important characteristics of the sys- 
tem. One of them is the structure of the intra-molecule 
Hamiltonian. At each molecule, we set an energy differ- 
ence between the states D{> 0) (see Fig.[T|) and different 
degeneracies of the states: gns and ^ls for the HS state 
and the LS state, respectively. We express the spin state 
at the i-th site by Si which takes —1 for LS and -1-1 for 
HS. The intra-molecule (on-site) interaction is expressed 



by 

Ho = IdJ2s^- (1) 

i 

If we take into account the effect of the degeneracy as 
a temperature dependent field, we can use an effective 
Hamiltonian with non-degenerate variables ai = ±1: 

n,s = lY.^D-~kBT\iig)a,, (2) 

i 

where g ~ gus I 5ls denotes the degeneracy ratio between 
the HS and LS states. 

The other important characteristic is the intermolecu- 
lar interaction. For the cooperative property in the SC 
transition, until recently a short-range Ising-type inter- 
action has been adopted in the so-called Wajnfiasz-Pick 
(WP) model:^ 

n = -jY,a,a,^\^Y.^D-kBT\x,g)u,. (3) 

This type of model has successfully explained various as- 
pects of the ordering processesi^Siiiii^ii^iiii However, the 
origin of the interactions between the spin states has re- 
mained unclear. There are various plausible origins of 
the interaction. 

As a possible interaction mechanism, the im- 
portance of elastic interactions has been pointed 
outpi^ii^'iiii^'i^'^aiSiiS^ The elastic constants may depend 
on the neighboring spin states. This dependence causes 
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an effective interaction between the spin states. This 
effect of the elastic constants was investigated in a one- 
dimensional (fD) two- level mode P^i^^ and also in a ID 
vibronic coupling model. In these one-dimensional ver- 
sions of the model, the clastic interactions can be traced 
out locally, leading to an exact mapping onto a ID 
Ising ferromagnet, so that there is no phase transition 
at nonzero temperaturesi22i2^ 

In higher spatial dimensions, as depicted in the inset 
in Fig. [1] the volume change of a molecule causes a dis- 
tortion of the lattice. Elastic interactions mediate the 
effect of this distortion over long distances. Therefore, 
in higher dimensions, the elastic interactions cause in- 
trinsically different effects than in one dimension. We 
denote this long-range interaction by 7iciastic({ci})- We 
do not know the explicit form of this interaction. (But see 
discussion in Appendix B). However, we recently demon- 
strated that this type of elastic interaction can induce 
a phase transition in spin-crossover systems^^-*^^'^- This 
elastic interaction model is a kind of compressible Ising 
model)^ and similar models have been studied for binary 
alloys»>^ 

Because the interaction originating from the elastic dis- 
tortions is qualitatively different from that of the nearest- 
neighbor Ising model, we are interested in the critical 
properties of systems with this type of interaction. We 
have previously studied phase transitions and the tem- 
perature dependence of ordering of model SC materials 
with specified parameters D and g. In those cases, most 
systems exhibit a first-order phase transition, and the 
critical properties of the models were not studied in de- 
tail. In the present study, we investigate properties near 
the critical point in the parameter space. In the case 
of the WP model, the critical properties are those of 
the short-range Ising ferromagnet. However, the criti- 
cal properties of the present model, i.e., the critical ex- 
ponents which characterize the critical universality, are 
expected to be different from those of the short-range 
Ising model. 

The organization of the rest of this paper is as fol- 
lows. In Sec. |TT]we present the model and the computa- 
tional method; in Sec. lIIII we discuss the finite-size scaling 
analysis of the critical properties; in Sec. IIVI we discuss 
the spin configurations and correlations; and in Sec. |V] 
we present a summary and discussion. A discussion of 
the long-range Husimi-Temperley model is given in Ap- 
pendix A, and a summary of finite-size scaling relations 
for mean- field phase transitions is given in Appendix B. 



II. MODEL AND METHOD 

In this paper, we study the critical phenomena of mod- 
els with elastically mediated spin-spin interactions on 
the simple square lattice (2D), and also on the sim- 
ple cubic lattice (3D) with periodic boundary condi- 
tions. Here we use Monte Carlo (MC) simulations ac- 
cording to the constant-pressure method. In the Monte 



Carlo simulation, we choose a site i randomly and up- 
date the spin state ai — (±1) and the position of the 
molecule(xi, j/i, Zi) by the standard Metropolis method. 
We repeat this update N times, where N is the number 
of lattice sites. Then, we update the volume of the total 
system. We define this sequence of procedures to be one 
Monte Carlo step (MCS). 

Instead of the Ising-like interactions of the WP model, 
Eq. ([3]), we adopt the following elastic interactions be- 
tween molecules)^ 

V = Kn + K„„ (4) 

Kn = YY.ln,-{R. + R,)r (5) 

Knn = y E [^*^- " ^(^'^ + ^j)]"' (6) 

where r^j is the distance between the i-th and j-th 
sites. Vnn expresses elastic interactions between nearest- 
neighbor pairs ((i,j)). Here, Ri and Rj are the radii of 
the molecules. The radius of each molecule is i?HS and 
i?Ls for the HS and LS states, respectively. In the present 
work, we set the ratio of the radii as Rhs/Rls — l-l- 
Vnnn cxprcsscs the elastic interaction of next-nearest- 
neighbor pairs (((i,j))), which is necessary to main- 
tain the lattice structure but not essential for the crit- 
ical behavior. We set the ratio of the elastic constants 
ki/k2 = 10. We set ki = 40 through out the present 
work. In this study, in order to exclude other effects 
than those due to elastic interactions through distortion, 
we assume that the stiffness constants ki and k2 do not 
depend on the spin state. If we were to allow spin depen- 
dence of ki and fc2, an effective short-range interaction 
would appear. In this sense, the present model treats 
only elastic interactions. 

The order parameter for the present model is the frac- 
tion of HS molecules, /hs = — !)• Hereafter, 
for convenience, we adopt the "magnetization," 

^ N 

i 

as the order parameter. In Fig. [21 we depict the temper- 
ature dependences of (M) for several values of D. Here, 
we find the typical Z?-dependences of {M{T)). That 
is, we find a smooth dependence for large values of Z3, 
and a first-order phase transition for small D. Between 
them, we have a second-order phase transition. This D- 
dependence is understood from the phase diagram of the 
non-degenerate model (i.e. g = 1)}^ In the present non- 
degenerate model a ferromagnetic phase transition takes 
place at T^^^ , and we expect a phase diagram as shown 
in the inset. In this phase diagram, H is the symmetry- 
breaking field. 

The temperature dependences of the state of the 
present model with degeneracy g > 1 (in this work we 
use g = 20) are given by the dotted lines in the phase 
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<M/N> 



FIG. 1; (Color online) Schematic picture of the energy struc- 
ture of a molecule. The left (right) minimum corresponds to 
the LS (HS) state. In the inset, schematic pictures of a lat- 
tice of LS molecules (left), and the distortion caused by a HS 
molecule in a lattice of LS molecules (right) are illustrated. 



diagram, 



H{T) = - {D-kBT\ng) , 



(8) 



When D is larger than — ksTclng, the temperature 
dependence of (M) is smooth, while it shows a first-order 
phase transition when D < Dc- If we consider a specific 
material, the parameters D and g are given, and the tem- 
perature dependence of the state is given by one of these 
dotted lines. In most cases, the ordering is either smooth 
or discontinuous, and the critical properties have there- 
fore not yet been seriously considered. 



III. CRITICAL PROPERTIES 



-1 



discontinuous 



critical 



H h 




H 1 



H 1 1 1 1 1 1 T 



1.0 



T • • 



discontinuous 
, critical 
..-smooth 



coexistence 



..■•VNBG 
•■■ >c 



FIG. 2: (Color online) The order parameter {M)/N vs tem- 
perature T for D = 2, 1.51, and 1. The temperature depen- 
dences are smooth, critical, and discontinuous (hysteresis), 
respectively. The inset shows a phase diagram of the model 
in the (T, H) plane, where T^^'^ is the critical temperature 
of the model Tidastic. The temperature changes in the model 
are given by the dotted lines in this phase diagram. When the 
dotted line crosses the coexistence line denoted by the bold 
line, the system undergoes a first-order phase transition. 



point. Here 



denotes the thermal average, i.e. 



TrM^e-^^ 
Tre-^« ■ 



(11) 



We study the critical properties of the elastically in- 
teracting model along the coexistence line given by T = 
D/\ng, i.e. H — Q?^ For the WP model, the critical 
properties are those of the Ising model. 

We next study the temperature (i.e., T = D /k-QThig) 
dependence of (M^). The spontaneous magnetization nis 
and the susceptibility per spin x s-i'^ obtained from the 
relation 



(M2) 
7V2 



(9) 



where N ^ L'^ \s the total number of spins, and 
^ NksT' 

which is the susceptibility per spin above the critical 



(10) 



In Fig. [31 we depict the temperature dependences of 
(M^) /iV^, and x~^- We find a clear linear dependence of 
(M2) /N^ below T ~ 0.20 for the 2D model, and T ~ 0.51 
for the 3D model. This linear dependence indicates that 
rTZg (xTc — T and thus the critical exponent (3 — 1/2. 

In Fig. [31 we also find that vanishes linearly at Tc, 
which indicates 7 = 1. This set of critical exponents 
agrees with those of the mean-field universality class. 
The size dependence of the inverse susceptibility in Fig. [31 
is rather large, but we found similar size dependences of 
(M^)/iV^, and x~^ in the long-range Husimi-Temperley 
model discussed in Appendix A. This indicates that the 
properties shown in Fig. [31 are inherent to models in the 
mean-field universality class f^i^^i^Ei^i^i^iMiS^i^^iSl 



4 




(a) 



0.4 



0.2- 




(b) 



FIG. 3: (Color online) Temperature dependences of (Af^) /N^ 
and the inverse susceptibility ~ NksT / {M'^) . (a) A two- 
dimensional system (square lattice): {M^)/N^ for the linear 
sizes L — 10, 20, 30, and 50 are plotted by closed symbols of 
circle, square, diamond, and triangle, respectively; and 
are plotted by the corresponding open symbols (multiplied by 
20 for improved visibility), (b) A three-dimensional system 
(simple cubic lattice): {M^)/N^ for L = 8, 12, 16, and 20 
are plotted by closed symbols of circle, square, diamond, and 
triangle, respectively; and are plotted by the correspond- 
ing open symbols (multiplied by 5 for improved visibility). 
The dotted straight lines show the expected behaviors for a 
mean-field phase transition in an infinite system. 
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FIG. 4: (Color online) Temperature dependence of the Binder 
cumulant. (a) A two-dimensional system (square lattice). 
Data for L — 10, 20, 30, and 50 are plotted as circle, square, 
diamond, and triangle, respectively, (b) A three-dimensional 
system (simple cubic lattice). Data for L = 8, 12, 16, and 20 
are plotted as circle, square, diamond, and triangle, respec- 
tively. 



d = 2 and 0.51 for d — 3. The value of U4 at the crossing 
is universal and independent of the spatial dimension. It 
is in excellent agreement with the theoretical resultj^Sii^ 



U4 = l- 



r*(i/4) 

247r2 



= 0.27---. 



(13) 



where F is the Gamma function. In Appendix A, we show 
the Binder plot for the long-range Husimi-Temperley 
model, which gives the same fixed-point value. 



We estimated the critical temperature by analysis of 
the Binder cumulant)^ 

Plotted for different system sizes, this quantity has a 
crossing at the critical point. It has been extensively 
studied for the mean-field universality classi^ We depict 
the Binder plot in Fig. [H The crossings are consistent 
with the values obtained from (]VP)/N^: ~ 0.20 for 



B. Finite-size scaling 

Finite-size scaling is one of the most useful methods 
to extract critical properties for infinite systems from 
numerical data for finite systems. However, special 
caution must be used when considering transitions in the 
mean- field universality class, which do not obey the hy- 
perscaling relation 2(3 -\- ^ = dv that relates the critical 
correlation-length exponent v with the spatial dimension- 
ality d for transitions with nonclassical exponents.*^. Es- 
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sentially, lengths are not well defined in systems with 
mean-field phase transitions, and the linear system size 
L is replaced by the number of sites N as the funda- 
mental finite-size scaling variable. A particularly clear 
example is the long-range Husimi-Temperley model dis- 
cussed in Appendix A, in which every spin interacts 
with every other with a strength proportional to 1/iV. 
The finite-size scaling variable that replaces the stan- 
dard tL^/'^ is W^/^^i^i This corresponds to an ejfec- 
tive exponent v* = 2/(i,— different from the value of 
I' = 1/2, obtained from the Gaussian approximation.— 
An effective exponent for the correlation function on the 
large scales that are relevant for finite-size scaling is ob- 
tained from v* by the standard exponent relation rj* — 
l—^jv* = (4— (i)/2. Thus one expects the scaling expres- 
sion (M2) = L'^+'^-^'M^itL^/"') = L3d/2_A/(2(^^d/2)^ 

where M.^ is a scaling function. A summary of the mech- 
anisms that lead to these results is given in Appendix B. 

In Fig. [5] we demonstrate that the Binder cumulants for 
different L collapse onto a single scaling function when 
plotted vs tL'^/'^ , and in Fig. [HI we plot the finite-size 
scaling functions for (M^). In both cases we find good 
data collapse, both for d = 2 and d = 3. These finite- 
size scaling relationships are also seen in the long-range 
Husimi-Temperley model discussed in Appendix A. 



C. Phenomenological scaling analysis 

In order to determine the critical temperature and the 
exponent rf — 2 — 7/1^*, the so-called phenomenological 
Monte Carlo rcnormalization plot is often useful.''^ That 
is, we plot 



U4 



DLOG 



ln((Af2)^/(Af2 



(14) 



\tl{L/L') 

as a function of T. The data for different sets of L and 
L' are expected to cross at a point which gives Tc and 
^/v* — 2 — rj* . In Fig. [71 we plot the temperature de- 
pendence of this quantity for two- and three-dimensional 
systems. We find a crossing in each figure at the position 
estimated by the values obtained in previous subsections: 
in the two-dimensional case, 



Tc ~ 0.20 and 77* ~ 1, 
and in the three-dimensional case, 

Tc ~ 0.51 and 77* ~ 0.5. 



(15) 



(16) 



We find a similar dependence in the Husimi-Temperley 
model given in Appendix A. 



IV. SPIN CONFIGURATION 

A. Spin correlation function 

Here, we concentrate on the characteristics of the spin 
correlation function. In Fig. [51 we depict the size and dis- 
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FIG. 5: (Color online) Finite-size scaling plots of the Binder 
cumulant U4 vs (T — Tc)L'''^^. (a) The two-dimensional sys- 
tem. Data for L = 10, 20, 30, and 50 are plotted as circle, 
square, diamond, and triangle, respectively, (b) The three- 
dimensional system. Data for L = 8, 12, 16, and 20 are 
plotted as circle, square, diamond, and triangle, respectively. 



tance dependences of the correlation functions for various 
values of T: T — 0.7/ ln((7), which is in the paramagnetic 
phase, T = 0.6/ln(.g) w T^, and T = 0.5/ln(.g) < T^. We 
plot the correlation function along the diagonal direction, 
i.e., Cd(r) = {(T{xo,yo)a{xo + r,yo + r))M We find un- 
usual spin correlation functions in the disordered phase. 
In short-range interaction models the correlation func- 
tion decays exponentially. In contrast, we here find the 
correlation to be nonzero and almost constant at long dis- 
tances in the disordered phase at T = 0.7/ \n{g). This ob- 
servation indicates that the spins arc strongly correlated, 
even at high temperatures. In the disordered phase, the 
susceptibility is an extensive quantity, and thus the total 
sum of the spin correlation function must be proportional 
to N: 



N N 



(17) 
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FIG. 6: (Color online) Finite-size scaling plots of {M^). (a) 
A two-dimensional system (square lattice). Data for L — 
10, 20, 30, and 50 are plotted as circle, square, diamond, and 
triangle, respectively, (b) A three-dimensional system (simple 
cubic lattice). Data for L = 8, 12, 16, and 20 are plotted as 
circle, square, diamond, and triangle, respectively. 



In order to satisfy this property, the constant value of the 
correlation function at long distances, cq, must depend 
on the system size as 



(18) 



This is in stark contrast to the result for Ising models 
with short-range interactions, cq ~ e~^/^^ with a corre- 
lation length ^ of order unity. 

At the critical point {T = 0.6/ln{g) ^ Tc), the size 
dependence of cq is given by 



Co oc 



1 



(19) 



This constant component at the critical point was 
pointed out by Luijten and Blotci^^ These observations 
are qualitatively different from those of the short-range 
Ising model. In the ordered state (T ~ Q.5/\n{g)), cq 
is independent of which corresponds to spontaneous 
order. 
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FIG. 7; (Color online) Phenomenological Monte Carlo renor- 
malization plots of {M^)/N, as defined in Eq. (fT4)) . (a) A 
two-dimensional system (square lattice). Data for {L,L') = 
(10,20), (10,30), (10,50), (20,30), (20,50), and (30,50) are 
plotted as cross, circle, square, diamond, triangle, and plus, 
respectively, (b) A three-dimensional system (simple cubic 
lattice). Data for (L,L') = (8,12), (8,16), (8,20), (12,16), 
(12,20), and (16,20) are plotted as cross, circle, square, dia- 
mond, triangle, and plus, respectively. 



B. Spin configuration in equilibrium 



Next, let us discuss the characteristics of the spin con- 
figurations in the model. In Fig.[9l we depict three snap- 
shots of spin configurations (a) at a high temperature, 
(b) near the critical point, and (c) at a low tempera- 
ture. We find that there are no large domain structures, 
even near the critical point. For comparison, we depict a 
configuration at the critical point of the two-dimensional 
nearest- neighbor Ising model (d). The difference is strik- 
ing. We also found the structure factor to be almost 
wave-number independent (not shown). From these ob- 
servations, we expect that usual critical behavior asso- 
ciated with two-phase coexistence will be suppressed in 
the present model. 
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FIG. 9; (Color online) Snapshots of equilibrium configura- 
tions of the two-dimensional model at (a) T — 1.2Tc, (b) 
T = Tc, and (c) T = 0.8Tc. (d) A snapshot of the nearest- 
neighbor Ising ferromagnet near the critical point, T = 2.3J. 



0.5- 



10 



J 



20 



30 



FIG. 8: (Color online) The size and distance dependences of 
the spin correlation function in the two-dimensional model 
at several temperatures, (a) T — 0.7/ln{g), which is in 
the paramagnetic phase, (b) T — 0.6/ln{g) ~ Tc- (c) 
T = 0.5/\n{g) < Tc in the ordered phase. In these figures, 
the system sizes are L = 10 (x), 20 (•), 30 (o), and 50 (*). 




FIG. 10; (Color online) Time dependence of the magnetiza- 
tion at the end of a hysteresis loop. 



C. Spin configuration at the end of the hysteresis 

loop 

We also studied the change of the configuration at the 
end of a hysteresis loop, i.e., near the (pseudo)spinodal 
that marks the limit of the metastable HS phase. For this 



purpose, we decreased the temperature gradually from 
the HS phase in the two-dimensional model with D — 0.4 
(< Dc — 0.6). The HS state remains as a metastable 
state beyond the coexistence curve. However at a certain 
point, it relaxes quickly to the LS state, marking the 
end of the hysteresis loop near the (pseudo)spinodal.''^ In 
Fig. 1101 we plot the time dependence of the magnetization 



8 





(e) (f) 

FIG. 11: (Color online) Time dependence ol the spin con- 
figuration at the end of a hysteresis loop, (a) Configuration 
just before the end point, at T = 0.1040. Here no change 
occurs after 40000 MCS. (b)At T = 0.1024. At this tem- 
perature, a sudden change of the magnetization occurs at 
22 000 MCS. (c) configuration at 23 000 MCS, (d) configu- 
ration at 24 000 MCS, (e) configuration at 25 000 MCS, (f) 
configuration at 26 000 MCS, 



as we decrease the temperature in steps by AT = 0.008 
every 40 000 MCS. In the figure, the temperature is kept 
fixed at T = 0.1040, 0.1032, 0.1024 and 0.1016. A rapid 
change of phase takes place at T = 0.1024. In Fig. [11] we 
show configurations during this rapid change (denoted by 
circles in Fig. fTU)) . 

In contrast to short-range interaction models, in 
which the phase change occurs through nucleation and 



FIG. 12: (Color onUne) A spin configuration at 200 MCS 
after a sudden change of T (T = lOTc to T = O.STc) along 
the symmetry line, D = T\ng. 

growth of compact critical droplets of the bulk equi- 
librium phase, the present system remains macro- 
scopically uniform during the whole transformation 
process. This is consistent with the accepted pic- 
ture of spinodal nucleation in systems with long-range 
interactions ,'*^''*^i"^^'^^i^°'^^ where the critical droplet is 
known to be extended and highly ramified with a density 
close to that of the metastable phase. It is thus extremely 
difficult to distinguish from the metastable background. 
Growth of this critical droplet occurs by a fiUing-in of 
its "interior," which is seen as the uniform change in the 
order parameter in Fig. 1101 



D. Spin configuration after quench into the 
low-temperature phase 

In short-range models with non-conserved order pa- 
rameter, the cluster size increases proportionally to the 
square-root of the elapsed time after a sudden quench 
from a disordered phase to a low-temperature phasei^ 
In contrast, the present model does not show such clus- 
tering configurations. In Fig. 1121 we show a typical con- 
figuration after quenching. Here we again find no large 
cluster growth, which indicates that there is no critical 
opalescence in the present model. 

These processes keeping uniformity can be understood 
in the following way. If a large domain exists, it causes 
a large distortion of the lattice, which is energetically 
unfavorable. Thus the system tends to be uniform on 
large length scales. This mechanism would be a char- 
acteristic of the present elastically induced mean-field 
phase transition. Beside the present SC system, there 
are various systems in which elastic interactions play an 
important role. For example, for the martensite transi- 
tion in metalsf^ the elastic interaction is important, and 
we expect similar critical behavior there. 

As mentioned previously, as far as a specific material is 
concerned, D is given, and the temperature dependence 
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of ordering is given by the dotted lines in Fig. 2. Thus, in 
most cases the phase transition is of first order. In such 
cases, the D dependence of the ordering studied in this 
paper is difficult to observe. However, the fact that the 
system is always uniform and no clustering occurs should 
be observable, even in a specific material. Moreover, by 
making use of the pressure dependence, ■22. we may also 
observe the critical properties and confirm the mean-field 
universality class. 



V. SUMMARY AND DISCUSSION 

We studied the critical properties of the elastically in- 
duced spin-crossover phase transition, finding it to be- 
long to the mean- field universality class. The temper- 
ature dependences of the long-range order and the sus- 
ceptibility were obtained in two- and three-dimensional 
models, and the corresponding critical exponents (3 and 
7 were found to be 1/2 and I, respectively, in agree- 
ment with the mean- field universality class. The size- 
and temperature-dependence of (M^) converged onto a 
scaling function. In the analysis of the finite-size scaling, 
we need critical exponents for the spin correlations, i.e., 77 
and V. We found that the effective values, 77* = (4 — d)/2 
and v* =2/d, are good for the scahng plots, as has been 
pointed out in various studies of the mean-field universal- 
ity class. We also found that the critical properties of our 
model agree well with the long-range interaction model 
(Husimi-Temperley model) , in which the spin correlation 
function is constant at large distances. 

We also studied characteristics of the spin configura- 
tions of the present model with effective long-range inter- 
actions. We found that the system does not show config- 
urations with large clusters, even following sudden tem- 
perature quenches, or at the edge of the hysteresis loop 
near the (pseudo)spinodal. Thus critical opalescence and 
conventional nucleation phenomena do not appear in the 
present model. In materials, it is difficult to change D 
or but the pressure dependence of these parameters'^ 
will enable them to be controlled, and we hope that the 
characteristic behaviors uncovered in this study will be 
found in real experiments in the future. 
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Appendix A: Finite size properties of the 
Husimi-Temperley model 

We study the finite-size dependence of the magnetiza- 
tion of the Husimi-Temperley model as a reference of the 
mean-field type behavior. The Hamiltonian is given by 




(20) 



Following the standard method, we obtain the partition 
function: 



Z = Tie-f^^ = e~ l^'=i 



dx 



27r 



x^ + N\n [2 cosh (a;y^2/3j7/V) ] 



If we estimate this integral by the saddle-point method, 
we obtain the mean-field free energy 



MF 



N 



--z^ + In 



2 cosh (zv^W 



(22) 



Here, we obtain the physical quantities of the model for 
finite values of N . The average of the square of the mag- 
netization is given by 

(M2) = = iV^ In Z 



dp J 



: tanh 



VNdz 



00 ^2tt 



2 +N In [2 cosh (zy'2/3j) ] 



where 



N 



(23) 



(24) 



The temperature dependence of (A/^) /N'^ corresponds to 
the square of the spontaneous magnetization Trig, 



(25) 
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FIG. 13: (Color online) Temperature dependence of 
{M^)/N'^ and the inverse susceptibility = TN/{M'^) for 
the Husimi-Temperley model. The circle, square, diamond, 
triangle and inverse-triangle denote A'^ = 100, 400, 900, 1600, 
and 2500, respectively. The open and closed symbols denote 
{M^)/N and X~^i respectively. 



FIG. 15: (Color online) The scaling plot of (A/^) for the 
Husimi-Temperley model. The circle, square, diamond, tri- 
angle and inverse-triangle denote N = 100, 400, 900, 1600, 
and 2500, respectively. 
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FIG. 14: (Color online) The Binder plot for the Husimi- 
Temperley model. The circle, square, diamond, triangle and 
inverse-triangle denote iV = 100, 400, 900, 1600 and 2500, re- 
spectively. 



The size- and temperature dependence of (M^) is found 
to converge in the standard finite size scaling plot as de- 
picted in Fig. 1151 The values of the correlation functions 
at large separations in Eq. (|18p and Eq. correspond 
to the above size dependences. 

These figures qualitatively agree well with those for the 
model of the elastic interaction mediated spin-crossover 
materials. 

In Fig. [16] we plot the phenomenological scaling plot of 
the present data 

for various sets of {L, L'). Here, we define L — N^^^. In 
general, if we use a definition L = N^^'^, DLOG becomes 
DL0Gx((i/2). 

Appendix B: Finite-size scaling of the mean- field 
model 



and k-QTN/ {NP) corresponds to the inverse susceptibility 
X^^ above the critical temperature. We plot the data in 
Fig. 1131 In the present model, the critical temperature is 
2J /k-Q- Hereafter we put J — 1 and k-Q — 1. 

The Binder plot of this model is depicted in Fig. [TH 
We find that the data for large N show a good crossing 
at Tc. 

The following size dependences are easily obtained: 

r 7V2 T<Tc, 
(M2) oc <^ iV3/2 T = T„ (26) 
\ N T>Tc. 



In this Appendix we summarize the finite-size scaling 
properties expected in a spatially extended system with 
mean-field behavior, which agree with those observed nu- 
merically in this paper. 

A d-dimensional cj)'^ lattice field theory with interaction 
range TZ can be defined by the Ginzburg-Landau Hamil- 
tonian in reciprocal space, 

k 

+ '/'ki0k2 0k3'/'-(ki+k2+k3) 

ki,k2,k3 
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FIG. 16: The phenomenological renormalization plot for the 
Husimi-Temperley model. Data for {L,L') = (10,20), (10,30), 
(10,50), (20,30), (20,50), and (30,50) are plotted by cross, 
circle, square, diamond, triangle, and plus, respectively. 



'k=0 



(28) 



where N = L"^ is the number of lattice points, t ^ {T — 
Tc)/Tc, and h is an applied magnetic field. For a = 2 and 
constant 7?,, this model has local interactions and upper 
critical dimension du = 4. For d > 4 it has classical 
mean- field critical exponents, for d = 4 it has mean- field 
exponents with logarithmic corrections, and for d < 4 
it has nontrivial critical exponents corresponding to the 
d-dimensional Ising universality classJ^ 

For d below four there are several ways the model can 
be modified to show mean- field critical behavior. One 
is to keep a = 2 fixed and let 7?, — > oo while using a 
scaling ansatz equivalent to a Ginzburg criterion^S? as is 
often done in studies of crossover scaling.— In this limit 
of infinitely weak, infinitely long-ranged interactions, the 
model reduces to the Husimi-Temperley model discussed 
in Appendix A. However, the method most relevant to 
elastic systems^^i^iSi^ is probably to increase the inter- 
action range by modifying <7i^ For < ct < d/2, this 
lowers the upper critical dimension to 



d„(cr) = 2(7 



(29) 



and leads to classical mean-field critical behavior for 
d > du{cr). (As for a = 2, classical exponents with loga- 
rithmic corrections are found for d — du{cr).) 

^From the terms corresponding to k = in Eq. (pS)) . 
one gets the standard mean-field critical exponents for a 
spatially uniform system, 

/3 = 1/2 (30) 

for the temperature dependence of the order parameter, 
^ oc \tf for t < 0, 



for its field dependence, cx \h\^/^ for t = 0, and 

7 - 7' = 1 (32) 

for the corresponding susceptibility, x — dcf)/ dh oc \t\^'^ ■ 
Spatial fluctuations are governed by the k'^ term. In 
the Gaussian approximation this yields 



= 1/cr 

for the correlation length, f oc Itl^"^, and 
7] = 2- a 



(33) 
(34) 



for the spin correlation function, c(r) oc 
exp[— r/^]r~'^'^~'^+''^42i However, renormalization of 
the "dangerous irrelevant variable" u that multiplies the 
fourth-order term in Eq. ([28]) '^^'^^ causes the fluctuations 
on large length scales comparable to the linear system 
size L instead to be governed by the cr-independent 
effective exponentsr^i^ 



and 



= 2/d 



2 -d/2 



(35) 



(36) 



Using these effective exponents in standard finite-size 
scaling relationS)^!^, one obtains the following scaling 
relation, 



(M^) oc L 



2d 



(37) 



(5 = 3 



(31) 



where the scaling relation j/iy* — 2 ~ rf has been used, 
and M'^{x) is a scaling function. ■^i The Binder cumulant 
Ui also becomes a scaling function of a; = tL'^/'^, ranging 
from 2/3 for a; ^ to for a; with the fixed-point 
value of Eq. (jl3p at a; = 0. Similarly, the phenomenolog- 
ical renormalization plot obtained from Eq. p4p will go 
from 2(ifor< < through d+ 2 -77* = 2{d~P/v*) = 3d/2 
at i = 0, to —d for t ^ 0. The spin correlation function 
at r - L takes the value cq oc ^-(^-2+')*) = L''^/'^. 
For r ^ L, on the other hand, the behavior is expected 
to be governed by the Gaussian exponent 77. However, 
much larger systems than the ones studied here would 
be needed to detect this behavior, which would enable 
one to measure the value oi 

We finally note two interesting aspects of these mean- 
field finite-size scaling relations. First, by making the 
replacement L'* = iV, it is easy to see that they all be- 
come independent of d, as long as d > du{cr). Second, we 
note that the Gaussian exponents can be used as in ordi- 
nary finite-size scaling theory, provided that L is replaced 
by the modified system size, 

The details of the effective long-range interactions in- 
troduced by the elastic degrees of freedom in the present 
system are not known, except for d = 1. In this case it 
has been shown rigorously that the model can be mapped 
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onto an Ising chain with nearest-neighbor ferromagnetic 
interactions and thus it exhibits no phase transition at 
nonzero temperatures. Much work has been devoted over 
the years to the interactions between defects in three- 
dimensional elastic solids, and it is generally argued that 
the dominant long-range interactions are dipole-dipole 
interactions ~ 1/r^ that can be attractive or repulsive 
depending on the relative orientations of the dipoles.^^ 
The case of d = 2 has been much less studied, and rele- 
vant works are much more recent. Dimensional analysis 
indicates that dipole-dipole interactions ~ 1/r^ should be 
present unless forbidden by symmetryi^iSii^ Although 
the effects of distortions in the present model arc not 
identical to those in the classic elastic media for which 
these results were obtained, we do not think it is unrea- 
sonable to assume that the elastically mediated interac- 



tions in our model are of such long range type. This then 
would lead to tr = and consequently d > ci„ = for 
both cases, so that classical mean-field critical behavior 
would indeed be expected. However, a weaker condition 
on cr, which still would lead to mean-field critical behav- 
ior, is obtained by simply requiring d > (i„(cr), leading 
to < cr < d/2 or, equivalently, interactions ~ i/j-d+c^ 
with d < d + a < 3d/ 2. We also note that the mecha- 
nism involving a = 2 and TZ oo leads to the same ef- 
fective exponents as the variable-tr mechanism,— 122 and 
so it would also be consistent with our numerical data. 
Which (if any) of these mechanisms best describes the 
long-range interactions that cause the mean-field critical 
behavior observed in the model studied here, remains an 
interesting question for future research. 
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